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Quantum walks provide simple models of various fundamental processes. It is pivotal to know 
when the dynamics underlying a walk lead to quantum advantages just by examining its statistics. 

A walk with many indistinguishable particles and measurements of non-classical multi-particle cor¬ 
relations is likely to reveal the quantum nature. The number of elements 0{n ) in a setup realizing 
walks grows with their length or spread n. We introduce the concept of a quantum leap, a process 
which can be achieved with fewer or complementary resources and which in a single step simulates 
another long process. The process and its leap are described by the same Hamiltonian but, the 
latter parametrizes the evolution with a tunable parameter of a setup. In the case of walks, a leap 
immediately gives a probability distribution which results only after many steps. This may be ap¬ 
pealing for simulation of processes which are lengthy or require dynamical control. We discuss a leap 
based on the multi-particle Hong-Ou-Mandel interference, an inherently quantum phenomenon. It 
reproduces a quantum walk enabling perfect state transfer through spin chains. It requires a beam 
splitter, two detectors and n particles to mimic a walk on a chain of size 0(n), for time fixed by 
beam-splitter’s reflectivity. Our results apply to a broad class of systems where the HOM-like effects 
can be observed, and may constitute a new approach to simulation of complex Hamiltonians with 
passive interferometers. 


Introduction. Quantum optics offers precise control 
and manipulation of low-energy quantum light and mat¬ 
ter. This is the first step towards efficient simulation 
of processes that are beyond the reach of present-day 
experiments and implementation of quantum-enhanced 
technologies. Here quantum walks (QWs) are receiv¬ 
ing much attention since they provide simple models of 
various fundamental processes in nature ranging from 
chaos [1], topological phases [2, 3] or photosynthesis [4- 
] to universal quantum computation [8] and quantum 
search algorithms [9-11] . Recently, numerous elementary 
QWs have been performed with photons [12—21] , trapped 
ions [22, 23], neutral atoms [24-26] and interacting atoms 
in optical lattices [27]. However, verification of the quan¬ 
tum nature of QWs is challenging. Also the number of 
elements in a setup realizing QWs grows with their length 
or spread. These facts are unfavorable especially for long 
walks, where decoherence plays a detrimental role. To 
circumvent these limitations, an alternative way to the 
future technological advances in QWs implementation is 
asking: are there genuinely quantum processes which ac¬ 
celerate QWs’ evolution, and which can be achieved with 
fewer or complementary resources? Here we examine one 
such example. 

QWs are the quantum analogy of the classical ran¬ 
dom walk in which an initially localized particle (a 
walker) evolves into a superposition over many positions 
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FIG. 1. Multiphoton Hong-Ou-Mandel interference imple¬ 
ments a quantum leap. The probability distribution of the 
photon population difference A in the exit modes of a beam 
splitter simulates a continuous time quantum walk obtained 
after longer evolution. 


in space [28-30]. Usually, the space is discrete and can 
be represented as a graph. The distinctive property of 
a QW is U-shaped statistics of the walker’s positions in 
a graph with variance scaling quadratically in time (bal¬ 
listic spread). This is in stark contrast to a classical 
random walk where the probability distribution is Gaus¬ 
sian. However, there exist both, classical processes which 
mimic the QW statistics (e.g. classical waves explain the 
QW of a particle on a line [31, 32]) and QWs character¬ 
ized by the binomial distribution (e.g. QW with inho¬ 
mogeneous time-dependent coin [33]). Thus, it is crucial 
to know when the dynamics underlying a walk is useful 
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for quantum tasks, just by examining its statistics. QWs 
involving many indistinguishable particles are likely to 
reveal quantum nature, possible to verify by non-classical 
multi-particle correlations. 

QWs with two particles were demonstrated with pho¬ 
tons [15, 20, 21, 34] and atoms in optical lattices [27], 
where generalized Hong-Ou-Mandel correlations were 
used for testing classical inequalities. Nonclassicality of a 
multi-mode two-photon qutrit was inferred from selected 
two-mode intensity correlations. Certification of multi¬ 
partite correlations solely based on one- and two-body 
ones is a difficult task, and is proven only for specific 
qubit cases [35]. Photonic implementation of QWs on a 
graph of size O(n) used arrays of 0{n 2 ) beam splitters 
(BSs) and 0{n ) detectors allowing t = n steps [20]. Al¬ 
ternatively, arrays of O(n) coupled waveguides and detec¬ 
tors were considered [15, 34] however, the current tech¬ 
nology allowed a detector behind every second waveg¬ 
uide. Optical lattices required 0(n) sites and in situ 
atom-resolved detection [27]. 

Here, by means of an example, we introduce a new 
concept of a quantum leap (QL), a process which in a 
single step simulates another long process. The process 
and its leap are described by the same Hamiltonian but, 
while the former describes the evolution in time, the lat¬ 
ter parametrizes it with a tunable parameter of a setup. 
For QWs, a QL immediately gives a probability distri¬ 
bution which results only after many steps. This may be 
appealing for simulation of processes which arc lengthy 
or require dynamical control. 

We discuss a QL employing the multi-particle Hong- 
Ou-Mandel interference [36] (HOM QL), shown in Fig. 1. 
This is an inherently quantum phenomenon, easily veri¬ 
fied by the HOM interferometric contrast. The HOM QL 
reproduces a continuous-time QW (CTQW) [29] enabling 
perfect state transfer through spin chains [37, 38] . It re¬ 
quires a BS of reflectivity r, two detectors and n particles 
to mimic a walk on a chain of size 0(n), for t fixed by 
r. It provides a new insight into the quantum-to-classical 
transition [14, 16] triggered by a partial distinguishability 
of particles. This architecture allows control over deco¬ 
herence during the evolution. 

Although our analytical results originate from 
quantum-photonic background, they apply to a broad 
class of systems and grant an intuitive access to many- 
body quantum dynamics. We analyze their feasibility 
and conclude with a discussion and open problems. 

HOM quantum leap. The HOM interference [36] is ob¬ 
served for two photons impinging on a BS described by 
the Hamiltonian 

H bs =a* l be i ‘ p + atfe- ilp . (1) 

The annihilation operators a and b correspond to the in¬ 
terfering modes and <p is the phase difference between 
the reflected and transmitted fields behind the BS. We 


assume <p = n. If the probabilities of reflection and trans¬ 
mission are equal, two identical photons will always leave 
together through the same exit port- photon bunching. 
Similar effects hold for multiphoton Fock-number states: 
photons go out only in certain configurations, e.g. such 
that the difference between the occupations of the exit 
ports is even, whereas an odd difference never occurs [39- 
4 ]. These interference effects, consisting of enhancing 
some events and canceling others, he at the heart of the 
QL generated by multiphoton HOM interference. 

A QW comprises at least one walker traversing a graph. 
A CTQW [29] takes place in the position Hilbert space 
where the basis states |Aj) correspond to vertices A,; in 
the graph. In the case of a QW on a line, the posi¬ 
tion can take any integer value. The evolution of the 
walker is generated by a Hamiltonian H which sets am¬ 
plitudes of hopping from vertex A^ to A j per unit time, 
Qi.j = (Aj\H\Ai). They are zero for disjoint vertices. 
The walker, initially in a state |A), after time t is in 
a state |A out ) = e~ lHt \A), where the variable A out 
is governed by a doubly-peaked probability distribution 
p(A out ) with variance scaling quadratically with time. 
Even in an infinitesimal time difference the walker can 
be transferred arbitrarily far from the starting position, 
albeit with low amplitude. 

We will show that for interfering Fock states \S — N)= 
1°) and \ N ) = ^7nC l°)> h bs generates a CTQW 

on a line of length S + 1 with specific jump amplitudes 
and its QL with length set by the BS reflectivity r. 

We parametrize the input \S — N) |jV) using the mode 
occupation difference A = S — 2 N and denote it |A). 
Behind the BS, the infinitesimal evolution turn |A) into 
the superposition 

Hbs\A) = 9a-2,a |A — 2) + ?a+2,a IA + 2). (2) 

This suggests the physical interpretation of A as the 
starting position of the walker in the graph, from where 
it jumps to the right |A + 2) or left |A — 2). Because the 
length of the jump is two in the position of the parti¬ 
cle, the evolution is confined to a region A e {—£, —S + 
2,..., S — 2, S}. The jump amplitudes 

9a,a-2 = 9a-2,a = 5\/ (S + A) (5 — A + 2) (3) 

indicate an important class of CTQWs. Hbs corresponds 
to the S x spin-matrix (in the Schwinger representa¬ 
tion) and the two-mode Fock input state can be associ¬ 
ated with a fictitious particle of total spin ^ and eigen¬ 
value of the S z component A. Under these mappings, 
the HOM interference describes rotation of the S z com¬ 
ponent, equivalent to a CTQW enabling perfect state 
transfer in linear spin chains of length S + 1, with inho¬ 
mogeneous couplings qA,A—2 [37]. 

The HOM QL results from parameterizing the pro¬ 
cess generated by Hbs oc S x with a BS reflectivity. The 
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FIG. 2. (Color online) Probability distributions of positions of a walker in the HOM QL performed with S = 50 photons, 
obtained for initial position of the walker a) A = 0, b) —30, c) —50, and reflectivity: r = 0.1 - green, r — 0.2 - red, r = 0.5 - 
blue, and r = 0.9 - gray. 


probability distribution of the final positions of a walker 
p(A out ) is governed by the unitary Ubs = exp {—iOHss}, 
where r = sin 2 9 controls length of the evolution. The 
spectrum of Hbs is harmonic and the evolution is pe¬ 
riodic. For r = 1 any input state |A) is mapped to the 
opposite vertex of the graph |—A). If r = 0 the unitary 
becomes the identity and a recurrence is observed: the 
QW returns to the initial state. For other values of r, 
the interfering Fock states are turned into a superposi¬ 
tion with population difference in the exit modes A out 
given by 


p(Aout) — A 



( 4 ) 


yAout j out | 

where A = + f a, f a, — : (r(l 


- r )) S (±==) 


f ^out 
J + 



k runs 


from max {0, A ° U *~ A } to min {/^, /| = 

(For the explicit form of the output superposition and 
an easy-to-read form of (4) see the Supplementary Mate¬ 
rial). Importantly, regardless the values of S and A, the 
distribution (4) is characterized by the ballistic spread 


Var(Aout) oc \ ^ + S') d 2 , (5) 


adjusted by the reflectivity r{0). (The proof is in the 
Supplementary Material). Thus, by fixing r, the multi¬ 
particle HOM interference reproduced the statistics of 
the discussed CTQW with the variance corresponding to 
an advanced stage of the evolution. This is the HOM 
QL. 

In Fig. 2 we plot (4) for various starting positions in 
the HOM QL performed with S = 50 photons. In Fig. 2a, 
the walker starts in the middle of the line, i.e. A = 0, 
and its walk spreads uniformly in both directions, thus 
the distribution develops symmetrically with respect to 
this point. This scenario takes place if the two interfer¬ 
ing Fock states are equal. If the initial position is shifted 
from the middle, A = —30 in Fig. 2b, the distribution is 
asymmetric for all r except for r = 1/2 and the walker 
tends to move in the initially preferred direction. This 
situation is simulated by two unequal Fock input states. 
The case when the walker starts in one of the most dis¬ 
tant positions in the graph, the ends of the line A = ±5, 


correspond exactly to perfect state transfer over a chain 
of interacting qubits as discussed in [37]. Fig. 2c plotted 
for A = —50, shows that the initial wavepacket prop¬ 
agates without distortion and reaches the other end of 
the line with unit fidelity. This is modeled by interfer¬ 
ence of a multiphoton Fock state with the vacuum state. 
The odd-even comb structure is revealed in all plots: the 
probabilities of measuring an odd difference are zero. 

Quantum-to-classical transition. In the presence of de¬ 
coherence, QW statistics turn to the classical one—the 
binomial distribution [42], We expect this effect also 
for the HOM QL, where the degree of coherence is con¬ 
trolled by the distinguishability of interfering particles. 
Initially all photons carried the same polarization. Now, 
one beam comes in a superposition of two orthogonal po¬ 
larization modes and b\_ 

—> cos y + sin y 6^. (6) 

y € (0,§) introduces weights between the modes, and 
‘tunes’ the distinguishability. The particles are fully 
indistinguishable if y= 0, whereas if y= : | they are maxi¬ 
mally distinguishable: |IV) and \S — N) are orthogonally 
polarized. Transformation (6) leads to the interference 
of | N) with a two-mode Fock state superposition 

( S n N ) 1/2 cosn y (sin y) s - N ~ n \n) \S — N — n) ± , 
instead of the single-mode Fock state |S — IV), as before. 

In Fig. 3 we show how the doubly-peaked distribution 
shown blue in Fig. 2a modifies if particles become distin¬ 
guishable. With increasing y, full cancellation of certain 
events is impossible and the two peaks gradually shift to 
the center of the line and merge to create the binomial 
distribution for y = Distinguishability mimics deco¬ 
herence because it leads to interference of the multipho- 
ton states with the vacuum state thus, it implements the 
usual model describing particle loss. It causes quantum 
particles to behave as classical ones, whose statistics can¬ 
not mimic quantum coherence and therefore we observe 
a transition from the quantum to the classical domain. 

Multi-dimensional HOM QL. Since a BS sees orthogo¬ 
nal polarizations independently, the above model gener¬ 
alizes to a two-dimensional case. This requires that the 
particle evolves independently in each direction, i.e there 
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FIG. 3. Modification of the probability distribution depicted 
in Fig. 2a in blue (A = 0, r = 1/2), for various degrees 
of distinguishability of particles: a) y = 7r/24, b) y = 7r/6, 
c) y = 7r/3 and d) y = 7r/2 (fully distinguishable). 


are no operations entangling the two polarizations. One 
simply considers two simultaneously interfering pairs of 
states: IS 1 — N) and | N) in horizontal and 15 — M) and 
| M) in vertical polarization. Thus now, each input beam 
is a two-mode beam, with fully distinguishable modes. 
This scenario differs from the previous decoherence model 
in that this time detectors are capable of distinguishing 
between orthogonally polarized photons. By considering 
orthogonal spectral modes instead of polarization, it is 
possible to extend the model to arbitrary dimension. 

Feasibility. Various experimental platforms approach 
regimes of parameters necessary for demonstration of our 
proposal. Atoms in optical lattices in the Mott insulator 
phase have fixed atom numbers per lattice site. In situ 
atom-resolved detection allows measurement of second- 
order correlations between the sites. These systems offer 
tunable beam-splitter-like interactions and direct scala¬ 
bility to larger particle numbers [27]. 

HOM QL was implemented with two atomic Bose- 
Einstein condensates (BEC) consisting of 10 4 atoms. 
They were prepared in a near-single-mode mixture of 
Fock states using spin dynamics and a tunable beam¬ 
splitter interaction was applied [43] . Measurements were 
performed with resolution of 20 atoms. U-shaped statis¬ 
tics without the parity structure was captured. Full 
demonstration would require single-mode input states 
and high-resolution detection. Purity of Fock states 
(losses) was not a limiting factor, since the HOM- 
interference contrast is immune to losses (exemplary dis¬ 
tributions for mixed Fock inputs and losses are shown in 
the Supplementary Material). 

Attempts to photonic demonstration employing bright 
near-single-mode squeezed vacuum twin beams (SV) are 
reported in [44], SV were produced in parametric down 
conversion, followed by mode filtering. Improvement 
of this experiment requires the creation of Fock states. 


They could be conditionally prepared from two-mode 
SV states, which posses perfect photon-number corre¬ 
lations between the modes (the signal and idler). The 
idler would be measured by a photon counting detector, 
whereas the signal would end up in a known Fock state. 
For the HOM QL, one would need two such sources and 
detectors. The range of S = 10 would match best su¬ 
perconducting transition-edge sensors (TESs) [45] with 
photon-counting efficiencies near 100% and well-resolved 
photon-number peaks up to around ten photons [46]. A 
BS of a controllable reflectivity can be implemented by 
a Mach-Zehnder interferometer with a dynamically con¬ 
trolled phase shift applied to one of the inputs [47, 48]. 
(Fig. 2 for S = 10 is shown in the Supplementary Mate¬ 
rial). 

Conclusions. Quantum leaps may be an important 
concept in the framework of quantum simulations. They 
quickly restore long evolution by mapping physical time 
to a controllable parameter of a setup while preserving 
an original Hamiltonian. This is appealing for simulation 
of processes requiring dynamical control, e.g. transport 
in presence of dynamical disorder, using passive setups. 

The HOM interference of many indistinguishable parti¬ 
cles provides a link between a walk of a quantum nature 
(useful for quantum tasks) and a QL. It is character¬ 
ized by the second-order visibility, which for two identi¬ 
cal Fock states impinging on a symmetric BS is quantum 
-y( 2 ) = 1/(2— 1/n) >1/2 and cannot be explained within 
classical theory (the case of two non-identical Fock states 
is discussed in the Supplementary Material). Moreover, 
it imparts an effective nonlinearity leading to the ex¬ 
change interaction between photons at a BS [49]. Up 
to date, nonlinear effects and interacting bosons have 
not been explored in photonic quantum walks. Impor¬ 
tance of generalized HOM effects for practical purposes 
has been recognized [50]. Combining the HOM inter¬ 
ference with a standard QW architecture will constitute 
setup enabling observation of large-scale interference, bo¬ 
son sampling [51, 52], multi-particle entanglement and 
non-linear phenomena [49]. 

The HOM QL reproduces the statistics of a CTQW al¬ 
lowing perfect state transfer over spin chains. Its variance 
is set by reflectivity of the HOM BS and reveals ballistic 
spread. It replicates quantum-to-classical transition trig¬ 
gered by partial distinguishability of particles and estab¬ 
lishes a link between the fundamental indistinguishability 
of quantum particles and the wavelike coherent nature of 
the walk. We have discussed one-dimensional HOM QL 
but by exploiting various degrees of freedom, the scheme 
may be generalized to higher dimensions. 

Implementations of QWs without interferometers are 
gaining a lot of interest now [53]. Our approach rep¬ 
resents a substantial saving in elements and detectors - 
we reduce the complexity of the system to having a sin¬ 
gle BS, two detectors and 0(n) particles to implement 
a CTQW with n positions. This allows control of deco- 
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herence during the evolution, by shifting it to the initial 
state preparation and detection scheme. 

We have discussed the feasibility of HOM QLs with 
optical lattices, BECs and optical systems. Our insights 
are of a general relevance and apply to a broad class of 
bosons: diffracting neutron [54] and X rays [55], surface 
plasmon polaritons [56], and fermions [57, 58] where the 
HOM-like effects were observed. 

Our results lead to open problems. Is it possible to 
find interferometers giving QWs with arbitrary jump am¬ 
plitudes? One can consider e.g. QWs with interacting 
BECs or atoms in optical lattices with larger particle 
numbers per lattice site, leading to generalization of the 
results presented in [43] and [27]. How does the interac¬ 
tion influence the dynamics of the walk and the observed 
probability distribution? What classes of Hamiltonians 
can be simulated with passive interferometers, perhaps 
including complex input states? In this case, a single 
step of a QL would be understood as a direct mapping 
of the input to the output state. The intermediate oper¬ 
ations implemented in such setup would not necessarily 
correspond to successive steps of evolution. 

We believe that this work will establish a new approach 
to quantum simulations and to bottom-up study of the 
dynamics of complex systems. 
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APPENDIX A: HOM QL PROBABILITY 
DISTRIBUTION p( A out ) 


We will now consider an unbalanced beam splitter with 
reflectivity r, which transforms the input annihilation op¬ 
erators a and b in the following way: a r = \/l — ra+y/rb , 
a t = —y/ra + \/l — rb. Note that r is the single-photon 
reflectivity because it gives the probability of reflection 
of a photon in the case if a single photon impinges on the 
beam splitter. The incoming Fock state | K) a \L) b , where 

| K) = |0) and | L) = |0), is turned into the 


following superposition 


IVw) = 7Wm.^ 1 ~ ra ^ ~ 'F' a t) K (Vf a i + Vl - ra\) L |0) 


K L 


V~r 


k—0 1=0 
K-k+l 


(-i) 


K-k 


n - L-l+k xk+l jK-k+L-l 

VI — r a' r a} |0). (7) 


We are interested in the probability distribution of mea¬ 
surement outcomes of photon counting detectors located 
behind the beam splitter \{p,q\il>out)\ 2 ■ Two Fock states 
| K) and | L) can always be parametrized by the total 
photon number S = K + L and the population difference 
A = K — L. In this case \K) = |^±A) and \L) = 

Since p + q = S as well, and p — q = A out the probability 
distribution reads 


p(Aout) — 


^ S— A ou t ^ j ^ ff+Aput ^ i 

i £ 

E 


A-A„ 


S -A 5+A out 
2 ' 2 


fc=max ( 0 


A 0 ut-A 


S-A 

2 

k 


(r(l-r)) d (l==)- 
) (^)‘ 


5+A 

2 

>S'~l~2\out 


( 8 ) 

The above distribution corresponds to Eq. (4) in the main 
text. It is shown for S = 10 in Fig. 4 below. 

In a similar way, one can compute the probabil¬ 
ity distribution for mixed Fock input states px = 
Ek =0 ( K k )v K ~ k 0-v) k \kM and p L = £f =0 ^)v L ~ l (l- 
rj) l \l)(l\, with purity equal Trp 2 , and lossy detectors. 
Fig. 5 displays the distributions from Fig. 4 computed for 
input Fock states with purity 0.83, whereas Fig. 6 shows 
the distributions from Fig. 2 in the main text computed 
for input Fock states with purity 0.47. Fig. 9 shows an 
array of the distributions obtained for S = 10 and various 
purities of input Fock states and amount of losses. 


APPENDIX B: BALLISTIC SPREAD OF p(A 0 ut) 

At first let us note that the beam splitter Hamiltonian 

H B s = o^be iip + (9) 

(Eq. (1) in the main text), for ip = tt, corresponds to the 
Schwinger representation of the S x spin--| matrix, de¬ 
fined as S x = ° tb + abt . Secondly, a two-mode Fock state 
IS — IV) |IV) can be associated with a fictitious particle 
of total spin j and eigenvalue of the S z component equal 
to f, where S z = ata ~ fetfe . Now |S-IV)|IV) = |A). 
Here A = S — 2 N denotes the initial mode occupa¬ 
tion difference. The evolution generated by the beam 
splitter’s Hamiltonian Hbs oc S x is parametrized by 
the beam-splitter’s reflectivity r, Ubs = exp{— idHss} 
where r = sin 2 9. 
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FIG. 4. (Color online) Probability distributions of positions of a walker in the HOM QL performed with S = 10 photons, 
obtained for initial position of the walker a) A = 0, b) —4, c) —10 and reflectivity: r = 0.1 - green, r = 0.2 - red, r = 0.5 - 
blue, and r = 0.9 - gray. 



FIG. 5. (Color online) Probability distributions from Fig. 4 computed for input Fock states of purity 0.83, S = 10 a) A = 0, 
b) -4, c) -10. 


Under these mappings, the HOM QL is in one-to-one 
correspondence with a spin--^ particle whose evolution is 
given by the Hamiltonian proportional to the S x opera¬ 
tor. In this picture the position operator (of the walker 
in the graph) is proportional to the S z operator and the 
states with well defined position correspond to the eigen¬ 
states of S z : !§,■§■)■ It is therefore straightforward to 
show that in the Heisenberg picture the S z evolves as 

S z (9) = S z cos# + S y sin#. (10) 

Note that S y is a linear combination of rising S+ = a^b 
and lowering S- = b^a operators, S y = l ( s -~ s +'> ; there¬ 
fore for the states |f, ^) the average value of S y is zero 
and we obtain 

(s z {0)) = {s z ) cose. (ii) 

Moreover, 

S 2 z {6) = S 2 Z cos 2 9 + S 2 y sin 2 9 + {S y , S z } cos 9 sin 9. (12) 

Due to the same reason as above, for the states |f, -f) 
the average value of S 2 {9) equals 

(S 2 (0)) = ( Sl) cos 2 9 + (S 2 ) sin 2 0. (13) 

In addition, since |f, ) are eigenstates of S z , we have 

( S 2 ) = ( S z ) 2 . The value ( S 2 ) can be evaluated in the 


following way. At first, we note that 

c2 _ -(S- - S+) 2 - S 2 _ - S 2 + S+S _ + S-S+ 

s y ~ 4 4 

(14) 

For the states |f, -f) we have 

{S 2 y) _ (S+S-) + (A_^ + ) _ f(f + 1) ~ (#) 2 ^ (15) 

Therefore 

Var[S z (0)] = ( S 2 (9))-(S z (9 )) 2 = sin 2 e . 

(16) 

The variance oscillates because the HOM QL takes place 
on a finite graph. However, for short times the bound¬ 
ary effects cannot play any significant role. Indeed, by 
approximating sin 9^9 we get 

Var[^(0)]^i(^^ + 5)0 2 , (17) 

which confirms ballistic spreading. In the main text we 
denote Var[S z (0)] by Var(A ou t)- Fig. 7 shows the vari¬ 
ance of the probability distributions of final positions of 
a walker in the HOM QL as a function of the time of the 
walk, for various values of S and A. 
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FIG. 6. (Color online) Probability distributions from Fie. 2 in the main text computed for input Fock states of purity 0.47, 
S = 50 a) A = 0, b) -30, c) -50. 




9 9 



e e 


FIG. 7. Variance of the probability distributions of final po¬ 
sitions of a walker in the HOM QL for a) & c) S = 50 and 
the initial position of the walker A = 0 green, A = —30 - 
red and A = —50 - blue b) & d) S = 10 and A = 0 - green, 
A = —4 - red and A = —10 - blue. 


APPENDIX C: SECOND-ORDER VISIBILITY IN 
THE HOM INTERFERENCE OF FOCK STATES 


Fig. 8 shows the range of m and n for which the visibility 
exceeds 1/2 thus, takes quantum values. 
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Visibility of the HOM interference at a beam splitter 
with reflectivity r and transmitivity t = 1 — r is given by 

= _ 2rt (- n a n b :) _ 

rt (: n 2 + n£ :} + ( r 2 + t 2 ) (: n a n b :) ’ 


FIG. 8. The range where the second-order visibility of the 
HOM interference v ® exceeds 1/2 for two non-identical Fock 
states | n) and |m), computed for n, m < 10 (left column) and 
n,m < 50 (right column) and a) r = 0.36, b) r = 0.43, c) 
r = 0.39, d) r = 0.45, e) & f) r = 0.5. 


where n a = a/a, n b = b^b, the annihilation operators a 
and b correspond to the interfering modes and :: denotes 
the normal ordering. If v^ > 1/2, the HOM interfer¬ 
ence does not have explanation within classical theory. 
Importantly, the visibility is immune to losses. Thus, its 
value is the same for two interfering pure and mixed Fock 
states and can be verified with lossy detectors. For two 
non-identical Fock states | n) and | m) it equals 

2 rt nm 

v ( 2 ) =-. ( 19 ) 

rt (n 2 — n + m 2 — m ) + (r 2 + t 2 ) nm 
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FIG. 9. (Color online) Probability distributions p(A ou t) from Fig. 4 (r = 0.1 - green, r = 0.2 - red, r = 0.5 - blue) computed 
for various purities of input Fock states and amount of losses. 5 = 10 a) A = 0, b) —4, c) —10. 
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